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In the present paper we calculate the dielectric functions of an ultra-relativistic plasma, such as an 
electron-positron or a quark-gluon plasma. We use classical transport theory and take into account 
collisions within the relaxation time approximation. From these dielectric functions we derive the 
• dispersion relations of longitudinal and transverse plasma waves. 
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I. INTRODUCTION 



In this paper, we will study the dielectric functions and dispersion relations for ultra-relativistic plasmas (mCT). 
The ultra-relativistic limit is relevant for the study of the quark-gluon plasma 0,0 or a hot electron-positron plasma, 
which is present in various astrophysical situations but might also be produced in the laboratory Q. Dispersion 
relations are of interest because they describe the propagation of collective plasma waves and provide information about 
■ physical quantities like screening lengths, oscillation frequencies, damping and particle production rates, transport 
fNJ | coefficients, and the equation of state Q- 

Dielectric functions and dispersion relations of ultra-relativistic plasmas have been calculated using transport theory 
and thermal field theory 0, Q . The new aspect of the present investigation is that we take into account collisions. 
For this purpose we start from the Boltzmann equation and use the relaxation time approximation for the collision term 
together with Maxwell's equations. We follow the non-relativistic derivation presented in Refs. [Tol[lH . Note that in a 
Oh! plasma, electromagnetic fields may be produced by external sources. These fields change the distribution and motion 
of the charged particles in the plasma, creating induced charges and currents, which also produce electromagnetic 
' fields. This system must be analysed self-consistently. 
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II. MAXWELL'S EQUATIONS 

To obtain dispersion relations we start from Maxwell's equations in the plasma (c = 1). For this purpose, we 
introduce the induced carge and current density (p,ji), the complex conductivity cry, and the dielectric tensor ey. 
Assuming a homogeneous medium these quantities are related to the electric field Ei and electric displacement Di in 
momentum space by 

ji(w,k) = (Xij(uj,k)Ej(Lj,k), (1) 

Di(cj, k) = ey(w, k)Ej(u, k). 
In addition we note that the conductivity and the dielectric tensor are related by 

-* % -* 

tij{oj, k) = Sij + —crij(oj, k). (2) 

Using these definitions, in the absence of external sources, Maxwell's equations in momentum space become 

(k x B)i = -w Eij (lo, k) Ej , (3) 
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k-B = 0, 

(k x E)i =uoBi, 

ki6ij(uj, k)Ej = 0. 

Eliminating B we obtain, 

- hkj - w%) Ej = 0. (4) 
We further specialize to an isotropic medium where the dielectric tensor can be decomposed: 

cy(w, fc) = P?e'(u;, fc) + P^(w, fe), (5) 

pT _ g fetjfr pL _ • 



Using J3J) we have, 



where 



ij l 3 £2 ' U £2 ' 

e*(w,£) = ip£e«, e l (u,k) = P±e ij . (6) 



(fc 2 -oA')J r = 0, (o;V)^ = 0, (7) 



(E T ) i = P^, (E^^P^Ej. (8) 

Non-trivial solutions are obtained from 

= 0, (9) 
wV(w,fc)-fc 2 = 0. 

III. KINETIC THEORY AND DIELECTRIC FUNCTIONS 

Let us consider an ultra-relativistic electron-positron plasma. We use kinetic theory to calculate the quantities e l 
and e in I0- The essence of kinetic theory is the assumption that the system can be described using a distribution 
function /(p, r, t) which gives the probability density to find a particle of momentum p at position r at time t. In 
order to obtain an expression for f(p, r, t) we must solve a kinetic equation. For the collisionlcss plasma it is common 
to use the Vlasov equation: 



df(p,r,t) -j dfjp, r, t) g 
dt dr 



E+-(vx B) 

c 



= 0. (10) 

op 



In order to solve this equation we make the approximation that the system is not far from equilibrium and expand 
around the equilibrium distribution function: 

mr,t) = f (p) + Smr,t), f (p)= cP/ t + 1 - (U) 

Note that we have neglected masses in the expression above by taking the ultra-relativistic limit (m <C T). The 
factor 4 in the equilibrium distribution fo(p) comes from the 2 spin states and from taking into account electrons and 
positrons in l|ll)fl at the same time, which is possible since the final results depend only on e 2 . In addition we define 
the particle number density 



N(r,t) = J dpf(p,f,t), N = J dpf (p), dp:=-^ 



d 3 p 
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In a weakly ionized plasma, we can include the effects of collisions by using a BGK collision term 

HHEt i-e ., 

we use a collision term on the right hand side of (|10|l of the form 



f(p,f,t)-N(f,t) 



hip) 



Jdp'foip') 



(13) 



where v is a velocity independent collision frequency. In perturbative QED the collision frequency between electrons, 
positrons, and photons (elastic 2 — > 2 scattering processes in the Boltzmann collision term) is of order e 4 (see e.g. (l2$). 
Since we are interested in investigating the effect of the collision frequency on the dispersion relations (see below), 
we will vary the parameter v from and 1 in units of m 7 , which is the only scale other than T in ultralrelativistic 
plasmas. (Note that e ~ 0.3 is not far from 1). 

The BGK collision term corresponds to an improvement of the relaxation time approximation for the collision term 
of the Boltzmann equation IS]. It simulates the effect of close binary collisions with substantial momentum transfer 
[TH ] and therefore it is particularly useful for describing collisions between ions and neutral gas molecules in a weakly 
ionized plasma. However, it can also be used as a rough guide for describing collisions between charged particles, such 
as in a hot electron-positron plasma. Note that the effect of long distance collisions with small momentum transfer 
can be considered by using a mean- field electric field E in 1(10(1 fl3T ] . 

We substitute (|ll(l into (|10fl and i|13|) and linearize in the deviation from equilibrium Sf to obtain, 



smr,t) = 



-ieE 



a/o(p) 

dp 



where 



If 

T) ■= — dp5f{p, r, t) , D(p) :=u)+iv — k- v{p) , v[p) = -. 
N J p 



(14) 



(15) 



Solving 1)14(1 and ((15(1 for Sf(p, r, t) we obtain an integral expression for the induced current: 



ji = e I dpvi{p)8f[p,r,t) 



ie / dpvi{p) 



p dp 



D-\p) 



No J - • " J p 

Using JQ| we extract an expression for (oj, k): 

cr^ = -ie 2 f dpvi(p)vj(p) d f° P D^ip) 



' ~iVo I dP /o(P )D ~ iP ' 



(16) 



dp 

+ Wo J dpv i {p)hip)D- x {p) J dp'v 3 {p>) d J0±D-\p') 



l-^-J dp"fo{p")D-\p"^ 



(17) 



Using J5J) and JHJ we obtain: 



e l = 1 + ^(kkjVij), 

t 3 i . . e l 

e = - H (an) . 

2 2uj y 11 ' 2 



(18) 



In the ultra-relativistic case we take v — 1 which means that in 1(17(1 the radial integral over p and the angular 
integrals decouple. Thus, in contrast to the non-relativistic case |10|, we are able to perform an analytic calculation 
of the dielectric functions. We define the integrals: 
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Ji 



dy 



!J 



7T 

eJ' + l 12' 

2 



X 



dx ■ 
i z — x 
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■h = 
J 3 = 



dx 

i z — X 

dx 

^ Z X 



-X(z), 
2-X(z) 



(19) 



where z = (u> + iv )/k and X(z) = 2 - zln (frrj- Using the definitions of f (p) and D(P) given in i|ll[) and (|15() we 
obtain: 



e ( = 1- 



e* = 1 



2e 2 T 2 
ir 2 Lok 

e 2 T 2 



h Ji + ttt 



7 2 

2 



2fc (i-tJ 3 )r 



ir 2 Lok 



Substituting in from l|19|) and defining the effective photon mass m 7 = eT/3 we obtain: 



e l {u,k) = 1 + 



e*(w,fc) = 1 



3m 2 



lo + iv to + iv + k 
1 ^ In 



3to 2 



1lo(lo + ii/) 



2k 
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u) -\-iv — k 
(uj + iv) 2 



k 2 



- 1 



ii/ lu + iv + k 

1 — — In - 

2fe w + iv — k 



lo + iv lo + iv + k 
1 — — m 



(20) 



(21) 



2A: 



lo + iv — k 



Let us discuss the analytic structure of these dielectric functions by comparing them with the collisionless results. 

(1) . In the collisionless case, v = 0, (|21|l reduces to the expressions found by using the Vlasov equation |(|, or the 
high-temperature limit of thermal QED 0,@- The correction to the vacuum value t l,t = 1 is proportional to the 
square of the effective photon mass, i.e., of order e 2 T 2 . 

(2) . Note that the longitudinal dielectric function does not follow from the collisionless case by simply replacing 
lo by lo + iv. As a consequence, the longitudinal part of the polarization tensor, II' = k 2 (l — e ) does not depend 
only on (lo + iv)/k, but also on iv/k separately. The transverse part of the polarization tensor, on the other hand, 
IT' = lo 2 (\ — e*), depends only on (lo + iv)/k. 



(3). Due to collisions, the dielectric functions (|21|l have explicit imaginary parts from terms containing iv. In 
addition, even for v = 0, there are imaginary parts coming from the logarithms, which are related to Landau damping. 

Apart from the constant factor m 7 , l|21() agrees with the results for the dielectric functions of a non-relativistic 
degenerate plasma as given in . This agreement occurs because the radial and angular integrals decouple for both 
the relativistic and non-relativistic degenerate plasmas, and the angular integrals, giving the functional dependence 
on lo and k, are identical in both cases. Note that for lo = 0, the longitudinal dielectric function becomes independent 
of the collision rate v. Therefore, the (square of the) Debye screening mass m 2 D = k 2 [e l (Lo = 0) — 1] = 3m 2 is also 
independent of v. 



IV. DISPERSION RELATIONS 



The equations © and (|21|l must be solved numerically. We scale all dimensionful variables by m 7 and produce 
numerical results for Lo/m-y versus k/rriy for various values of v/m 1 . Figs. (1-4) show the real and imaginary parts 
of the dispersion relations for transverse and longitudinal modes. We correctly reproduce the known results from the 
Vlasov equation (|10fl in the limit v = . These results can also be derived from quantum field theory, using lowest 
order perturbation theory in the high-temperature or hard thermal loop approximation Isl Il4l| . The real part of the 
dispersion relation is remarkably insensitive to v. The imaginary part, which is identically zero when v = 0, acquires 
a non-zero contribution for v =/= 0. 
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FIG. 1: Real part of the transverse dispersion relation (the parameter v runs from to 0.8 in steps of 0.2 with the smallest 
value being at the top of the graph) 
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FIG. 2: Real part of the longitudinal dispersion relation (the parameter v runs from to 0.8 in steps of 0.2 with the smallest 
value being at the top of the graph) 
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FIG. 3: Imaginary part of the transverse dispersion relation (the parameter v runs from to 0.9 in steps of 0.05 with the 
smallest value being at the top of the graph) 
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FIG. 4: Imaginary part of the longitudinal dispersion relation (the parameter v runs from to 0.9 in steps of 0.05 with the 
smallest value being at the top of the graph) 
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V. DISCUSSIONS AND CONCLUSIONS 

In the present investigation we have calculated the dielectric functions and the dispersion relations of electromagnetic 
plasma waves in an ultra-relativistic electron-positron plasma, taking collisions into account. Adding the BGK collision 
term to the Vlasov equation, we have derived analytic expressions for the longitudinal and transverse dielectric 
functions in the case of a small deviation from equilibrium. Using Maxwell's equations, the dispersion relations of 
longitudinal and transverse plasma waves follow numerically from the dielectric functions. Surprisingly, the deviations 
from the collisionless case are rather small, resulting only in a small shift of the dispersion relations towards lower 
energies. For example, the plasma frequency, given by = 0) = LOxik = 0), is reduced from m 7 at v — to 

only about 0.9 m 1 at v — m 7 15]. However, when collisions are included, the longitudinal dispersion at v ^ 
intersects the light cone oj = k, in contrast to the case of collisionless dispersion, where ujL,r(k) > k for all k. In the 
collisionless case the only damping mechanism is Landau damping, which occurs for uj 2 < k 2 , therefore plasma waves 
are undamped for v = 0. For v ^ the plasma waves are damped by the collisions, which produces damping rates 
7l.t(&) = -Imwi^ffc). The longitudinal wave is also Landau damped for momenta that satisfy Wi(fc) < k. 

It should be noted that we have introduced the collision rate v phenomenologically as a constant free parameter. In 
hydrodynamical models one normally assumes that local equilibrium is achieved after some momentum independent 
relaxation time which can be identified with the inverse of the collision rate. In principle, the collision rate could 
be determined pcrturbatively leading to a momentum dependent rate proportional to e 4 |l2j. However, it would be 
much more difficult to work with a momentum dependent collision rate. 

The effect of collisions in a quark-gluon plasma, e.g., on the damping of plasma oscillations and energetic quarks 
and gluons, has been discussed extensively in the literature using perturbative QCD as well as transport theory 
(see e.g. ^2 Our results for the dielectric functions as well as the dispersion relations are valid for chromo- 

electromagnetic plasma waves in a quark-gluon plasma, if we replace the effective photon mass w 7 by an effective 
gluon mass m g = gT^J(l + ?i//6)/3, where n/ is the number of quark flavors with masses m <^T 5J. However, this 
result is not fully consistent: neglecting the non-abelian terms of order g in the Yang- Mills equations contradicts the 
asumption that we can use a phenomenological collision term to model all collision processes, including higher order 
processes in g. At best, our phenomenological approach produces a rough estimate of the effect of collisions on plasma 
waves in a quark-gluon plasma. 
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